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We consider a spatial model related to bond percolation for the spread of a disease that includes 
variation in the susceptibility to infection. We work on a lattice with random bond strengths and 
show that with strong disorder, i.e. a wide range of variation of susceptibility, patchiness in the 
spread of the epidemic is very likely, and the criterion for epidemic outbreak depends strongly on the 
disorder. These results are qualitatively different from those of standard models in epidemiology, 
but correspond to real effects. 



The most commonly used models for the spread of an 
epidemic assume perfect mixing: i.e., all individuals are 
able to infect all others. In an inhomogeneous popu- 
lation the susceptibility to infection is replaced by its 
average. To account for spatial effects many workers use 
reaction-diffusion equations to describe traveling waves of 
infection . In both cases important phenomena are not 
reproduced by the models. For example, it is common ex- 
perience that in plant diseases |^,|| islands of susceptible 
individuals can be protected by a band of immune ones, a 
'firewall'. A map of the infected areas looks patchy. Even 
in human diseases such as AIDS these effects seem to be 
important Q . Such effects do not occur in perfect mixing 
or reaction-diffusion theories in which susceptibilities are 
replaced by averages. In this work we give a simple model 
for these phenomena using ideas taken from percolation 
theory. We show that for inhomogenous populations, the 
variability of susceptiblity to infection among individuals 
gives rise to qualitatively new effects. 

Percolation theory [^,^ is a natural way to think 
about populations which are not well mixed. Consider 
a spatially random distribution of susceptible individu- 
als which occur with probability p. A (site) percolation 
model for an epidemic is a lattice algorithm where, with 
probability p, a site is occupied and susceptible sites in- 
fect their susceptible neighbors. For diseases where some 
agents travel long distances, 'small-world' theory ex- 
tends the percolation model. In this approach a few sites 
are considered to be neighbors of distant sites. 

Percolation accounts for the existence of firewalls and 
islands because near the percolation threshold only part 
of the lattice belongs to the spanning cluster. However, a 
close look makes such an application suspect. Islands and 
firewalls occur only in a rather narrow transition region 
5p near the critical probability, pc- In nature it would 
be unlikely that a fine-tuning of parameters to be near 
Pc would occur. This kind of difficulty always plagues 
attempts to apply theories of critical phenomena in the 
natural world. Here, we show that introducing disorder 



in the susceptibility to infection solves this problem in an 
elegant way because it broadens the transition region for 
the outbreak, so that a generic epidemic could produce 
islands and firewalls. 

There is a further unexpected result of this study. By 
definition, an epidemic starts when a sick individual in- 
fects more than one other before she recovers, i.e. i?o > 1, 
where Ro, the reproduction number, is the mean number 
of infections produced by a site. The usual expression for 
Ro is xS{Q)t where x is the probability per unit time to 
infect neigbors, 5(0) the number which can be infected, 
and T the time to recovery. We will show that disorder 
can change this formula by a large amount. 

Our model is of the SIR (Susceptible, Infected, Re- 
covered) type. We use bond percolation: each site of a 
square lattice is occupied by an individual which can be 
infected by neighbors connected by bonds (four in our 
case) . To account for variability the probability per unit 
time of infection along a given bond is chosen from a dis- 
tribution f{x). After exactly r time steps any infected 
site recovers. We start with one / in the middle of the 
lattice and the rest S. We say that we have an epidemic 
when the / and R sites span a large lattice, i.e. reach the 
edges. A snapshot is shown in Figure 1. If the recovery 
time, T, is too small the epidemic will die out, and if r is 
large enough it will persist. A critical value of r plays the 
role of Pc- In Figure 2 we show the spanning probabil- 
ity, M , for a 256 x 256 lattice as a function of t for two 
different choices of f{x). M is the probability to span 
starting from a single site. We should note that M is 
the same as the mass of the infinite cluster in percolation 
theory i.e., the fraction of sites which belong to the 
spanning cluster. This is true because the probability to 
span starting at any given site is the probability that the 
site is a member of the infinite cluster. 

We have studied two classes of /'s which we call weak 
and strong disorder. An example of weak disorder is: 

Other functions with a rather narrow range of x give 
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similar effects, and the results do not differ very much 
from the case with no disorder ||]. 

Strong disorder refers to broad distribution such as: 

fs{x) ^C/X Xmtn < X < 1; C = 1/| h-l{Xmin)\ (2) 

Any strongly skewed distribution should give similar re- 
sults to those quoted here. Distribution functions of this 
type were studied [Q in random resistor networks, and 
are known, in that context, to give different behavior 
from the ordered case. 

The black portion of Figure 1 corresponds to sites that 
are protected by the weak bonds at the edge of the in- 
fected cluster, the firewall. Effects of this type occur only 
in the transition region of Figure 2. The transition region 
is very broad for strong disorder. Firewalls and large is- 
lands of uninfected sites arise purely from the statistics 
of the problem, as in any percolation problem. 

Ordinary bond percolation has some bonds present and 
some absent, which appears to be quite different from our 
approach. If x takes on a single value, the relationship 
was pointed out by Grassberger ||]: the probability to 
infect a neighbor before recovery is p = 1 — (1 — a;)"^. Thus 
p is the fraction of bonds completed in a given epidemic. 
If we have disorder we write: 

p=l- J{1- xYf{x)dx. (3) 

For a square lattice pc = 1/2, so we expect that when p in 
Equation || is greater than 1/2 we will have an epidemic. 
We verify this in the inset to Figure 2 where plot M 
against p for various choices of f{x) and compare them 
to ordinary bond percolation on the same lattice. 

We can discuss the threshold for outbreak of an epi- 
demic by using the usual SIR equations for the perfect 
mixing case: 

dS/dt = -xSI, dl/dt = xSI - I/t (4) 

Here there is a rate of recovery 1 /r rather than a fixed 
recovery time. However this is qualitatively the same. 
Clearly, we have an outbreak if Ro = xS{0)t > 1. 

If we convert Equation ^ into a spatial model by adding 
diffusion terms then the criterion for outbreak is of the 
same form H], where S'(O) now refers to the initial number 
of susceptibles which can be infected by a single infected 
site, a number of order unity. For the Grassberger model 
(or weak disorder) we have a similar looking result: we 
put p = 1 — (1 — xY" = Pc ~ 1/2. For small x, Ro = 

XT, = 0(1). 

Strong disorder is quite different. We can estimate the 
integral in Eq. (||) using the fs {x) of Eq. (^) for Xmin < < 
1. We get 

J {1 - xY fs{x)dx « [-F - \n[TXrmn)]/\ \n{Xmi7i)\ (5) 



where F is Euler's constant. Using this estimate we find 
Tc « 0.56/^Xmin- To estimate Ro we replace x by its av- 
erage, X. The naive estimate for the reproduction number 
at the threshold for outbreak gives: 

Rq ^ TqX OC l/[\/:2^mm 

1 ln(a;„ii„)|]. (6) 

This estimate for Ro can be arbitrarily large. The effec- 
tive Ro for a disordered lattice is very much smaller than 

TcX. 

This statement may seem reminiscent of the well- 
known fact Q that for diseases such as AIDS the naive 
formula for Ro does not work if there are very active 
agents, because active agents have an outsize effect in 
spreading the disease. However, the results are, in fact, 
in opposite directions. Here, the spatial correlations be- 
tween strong bonds make them less effective than they 
might seem to be. A clump of very infective bonds does 
not do as much damage as we might suspect since the 
infected sites are sharing the same victims. In the tran- 
sition region there are always bottlenecks to the spread 
of the infections, and thus random spatial correlations 
are always important. 

We can understand the broadening of the transition 
region due to disorder in the same way. Using the ex- 
pressions in Eq. (^) and Eq. (||) we find 5t/tc — 
\\n.{xmin)\5p- We can define the transition region, 5p, 
as the range over which M increases from, say, 0.1 to 
0.9; it is of order 0.05. We see that 5t/tc is large as 

Xmin ^ 0. 

In epidemics of human diseases some (usually a small 
number) of individuals travel long distances. We use 
a small-world lattice to account for these contacts by 
adding a small fraction, 0, of 'long' bonds connecting 
randomly chosen sites The result is shown in Figure 
3. Even a very small breaks up the large uninfected 
regions. The percolation threshold is lowered, and below 
the old Pc there are many small regions of epidemic con- 
nected by long bonds. However, the effects of disorder 
remain. Long distance contacts are not the same as per- 
fect mixing: firewalls still occur over a substantial range 
of parameters. 

To understand these results we extend the work of 
Newman and Watts to our case. The infected islands 
of Figure 3 can be considered to be nodes of a random 
graph whose edges are the long bonds. The graph per- 
colates if there are two islands per long bond ||ll[]; this 
occurs for some p < p,. 

We need to count the nodes of the random graph. 
To do this we use the average island size defined as 
fi = J2 s^«s/ J2 where ris is the number of clusters of 
size s. This is the correct average because it counts the 
probability that a given site (the end of a long bond) be- 
longs to a cluster [|j . Scaling theory gives n = K \p—pc \ 
where K is a constant, and 7 is a critical exponent equal 
to 43/18 = 2.39 for our case. The number of islands is 
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Ni = N/n, where A'' is the total number of sites. The 
number of long bonds that connect two sites that are 
parts of islands is B = N4>[l-{l-p)'^]'^ . Setting 2N, = B 
gives a criterion for percolation on a small world lattice 
with one adjustable parameter, K: 

cb^K{p,~pr{l-{l-pr)-^ (7) 

In Figure ^ we show the depression of the percolation 
threshold as a function of 0. We exhibit data as a func- 
tion of p rather than r. They are related by Eq. (|^). 
The threshold is obtained by requiring that an epidemic 
infects 20% of the lattice sites. (In this case, with long 
bonds, the notion of spanning loses its meaning). We 
also show a fit to the formula in Eq. (Q) with 7 and K 
as adjustable parameters. The best fit gives 7 = 2.40, 
remarkably close to the theoretical value quoted above. 
However, as we see in Figure (||) there are still islands 
is surrounded by a firewall, and all of the above consid- 
erations carry over in the strongly disordered case. In 
particular, the size of the transition region in terms of r 
is expanded 

Our model for percolation with strong disorder illus- 
trates two effects: there is a mechanism which produces 
patches of uninfected but susceptible individuals without 
fine tuning of parameters. And for this case Ro depends 
on the disorder. Whether strong disorder is a valid de- 
scription of nature needs to be determined by consider- 
ation of real epidemics. It is known, for example, that 
variability in susceptibility increases when epidemics re- 
cur. Our model may apply best to that case. 
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FIG. 1. Snapshot of the final cluster of recovered after 
the epidemic has died out on a 256x256 lattice with periodic 
boundary conditions. The cluster of recovered is gray, and 
the uninfected are black. A strong disorder {xmin = e"'^") 
bond distribution (g|) is used. 
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FIG. 2. The probability, M, that an epidemic started from 
a point spans the lattice as a function of infection time t for 
two different bond distributions. Strong disorder broadens 
the transition region in r. The transition region (e.g. the 
range 0.1 < M < 0.9) is about 12 times larger for strong 
disorder than for weak. For we use Xmax ~ 0.003, fa 
uses Xmin = 6^^^. Averages are over 2000 simulations on a 
200x200 lattice. [Inset] Data collapse of A4 as a function of 
p, compared with bond percolation. The bond distributions 
used are: (o) bond percolation, (+) fs with Xmin = e~^^, (x) 
fs with Xmin = e"^, (*) with Xmax = .001, and (A) /„ 

with Xmax = .01. 
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FIG. 4. Dependence of the percolation threshold on (o) 
compared with the scaling theory prediction on a 200x200 
lattice. 
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